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^ , The separation method developed earher by us [Nucl. Phys. A598 390 (1996)] to 

calculate and analyze solutions of the BCS gap equation for ^Sq pairing is extended 
and applied to ^P2-'^F2 pairing in pure neutron matter. The pairing matrix elements 
are written as a separable part plus a remainder that vanishes when either momen- 
tum variable is on the Fermi surface. This decomposition effects a separation of 



in 
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\^ ' the problem of determining the dependence of the gap components in a spin-angle 

<^ ■ representation on the magnitude of the momentum (described by a set of functions 

independent of magnetic quantum number) from the problem of determining the 
' dependence of the gap on angle or magnetic projection. The former problem is 

' solved through a set of nonsingular, quasilinear integral equations, providing inputs 

for solution of the latter problem through a coupled system of algebraic equations 
for a set of numerical coefficients. An incisive criterion is given for finding the up- 
per critical density for closure of the triplet gap. The separation method and its 
^ , development for triplet pairing exploit the existence of a small parameter, given by 

^ I a gap-amplitude measure divided by the Fermi energy. The revised BCS equations 

admit analysis revealing universal properties of the full set of solutions for pair- 
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ing in the absence of tensor coupling, referring especially to the energy degeneracy 
and energetic order of these solutions. The angle-average approximation introduced 
' by Baldo et al. is illuminated in terms of the separation-transformed BCS problem 

and the small parameter expansion. Numerical calculations of ^P2 pairing parame- 
ters and gap functions, with and without coupling to the ^F2 state, are carried out 
for pairing matrix elements supplied by (vacuum) two-neutron interactions that fit 
nucleon-nucleon scattering data. It is emphasized that ab initio evaluation of the 
in-medium particle-particle interaction and associated single-particle energies will 
be required if a reliable quantitative description of nucleonic superfluids is to be 
achieved. 



1 Introduction 



The study of superfluidity in infinitely extended nuclear systems has a long 
history [1-82], predating the 1967 discovery of pulsars [83], which were soon 
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identified as rapidly rotating magnetic neutron stars [84]. Interest in nucle- 
onic pairing has intensified in recent years, owing primarily to experimental 
developments on two different fronts. In the field of astrophysics, a series of X- 
ray sateUites (including Einstein, EXOSAT, ROSAT, and ASCA) has brought 
a flow of data on thermal emission from neutron stars, comprising both up- 
per limits and actual flux measurements. The recent launching of the Chandra 
X-ray observatory provides further impetus for more incisive theoretical inves- 
tigations. Realistic ah initio prediction of the microscopic physics of nucleonic 
superfluid components in the interiors of neutron stars is crucial to a quan- 
titative understanding of neutrino cooling mechanisms [85-95] that operate 
immediately after their birth in supernova events, as well as the magnetic 
properties, vortex structure, rotational dynamics, and pulse timing irregular- 
ities [96-99] of these superdense stellar objects. In particular, when nucleonic 
species enter a superfluid state in one or another region of the star, suppres- 
sion factors of the form exp(— A^/fcsT") are introduced into the expression 
for the emissivity, A^? being an appropriate average measure of the energy 
gap at the Fermi surface. On the terrestrial front, the expanding capabilities 
of radioactive-beam and heavy-ion facilities have stimulated a concerted ex- 
ploration of exotic nuclei and nuclei far from stability, with special focus on 
neutron-rich species [100,101]. Pairing plays a prominent role in modeling the 
structure and behavior of these new nuclei. 

In this article wc will be principally concerned with nucleonic pairing in the 
astrophysical setting. To a first approximation, a neutron star is described 
as a neutral system of nucleons (and possibly heavier baryons) and electrons 
(and possibly muons) in beta equilibrium at zero temperature, with a central 
density several times the saturation density po of symmetrical nuclear mat- 
ter [102,103,98,104,99]. The gross structure of the star (mass, radius, pressure 
and density profiles) is determined by the Tolman-Oppenheimer-Volkov gen- 
eral relativistic equation of hydrostatic equilibrium, consistently with the con- 
tinuity equation and the equation of state (which embodies the microscopic 
physics of the system). The star contains (i) an outer crust made up of bare 
nuclei arranged in a lattice interpenetrated by relativistic electrons, (ii) an in- 
ner crust where a similar Coulomb lattice of neutron-rich nuclei is embedded 
in Fermi seas of relativistic electrons and neutrons, (iii) a quantum fluid inte- 
rior of coexisting neutron, proton, and electron fluids, and flnally (iv) a core 
region of uncertain constitution and phase (but possibly containing hyperons, 
a pion or kaon condensate and/or quark matter). 

It is generally accepted [55,57] that the neutrons in the background fluid of the 
inner crust of a neutron star will pair-condense in the ^Sq state. Qualitatively, 
this phenomenon can be understood as follows. At the relatively large average 
particle spacing at the "low" densities involved in this region (p ~ po/10), 
the delocalized neutrons experience mainly the attractive component the ^Sq 
interaction. However, the pairing effect is quenched at higher densities, ~ po 
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and beyond, due to the strong short-range component of this interaction. 
By similar reasoning, one expects ^Sq proton pairing to occur in the quantum 
fluid interior, in a density regime where the proton contaminant (necessary for 
charge balance and chemical equilibrium) reaches a partial density Pp ~ Po/lO- 

The energy dependence of the nucleon-nucleon (NN) phase shifts in different 
partial waves offers some guidance in judging what nucleonic pair-condensed 
states are possible or likely in different regions of a neutron star. A rough cor- 
respondence between baryon density and NN bombardment energies can be 
established through the Fermi momenta assigned to the nucleonic components 
of neutron-star matter [57], which provide a crude measure of the maximum 
relative momenta attained in the nuclear medium. In (pure) neutron matter, 
only T = 1 partial waves are allowed by the Pauli principle. Moreover, one 
need only consider partial waves with L < 4 in the range of baryon density - 
optimistically, p < (3 — 4)po - where a nucleonic model of neutron-star ma- 
terial is tenable. The ^Sq phase shift is positive at low energy (indicating an 
attractive in-medium force) but turns negative (repulsive) at around 250 MeV 
lab energy. Thus, unless the in-medium pairing force is dramatically different 
from its vacuum counterpart, the situation already suggested above should 
prevail: S-wave pairs should form at low densities but should be inhibited 
from forming when the density approaches that of ordinary nuclear matter. 

The next lowest T = 1 partial waves are the three triplet P waves ^Pj, with 
J = 0, 1, 2, which are actually coupled by the tensor force to the triplet F waves 
with the same total angular momentum J. Of these channels, the first two are 
not good candidates for pairing. In the ^Pq state, the phase shift is only mildly 
attractive at low energy, turning repulsive at a lab energy of 200 McV, while 
the ^Pi phase shift is repulsive at all energies. On the other hand, the '^P2 phase 
shift is negative at all energies, and indeed is the most attractive T = 1 phase 
shift at energies above about 160 MeV. This exceptional attractive strength 
in the ^P2 state is due mainly to the spin-orbit interaction and in part to the 
coupling to ^F2 wave via the tensor component of the A^A^ force. A substantial 
pairing effect in the ^P2-^F2 channel may be expected at densities somewhat 
in excess of po, again assuming that the relevant in- vacuum interaction is not 
greatly altered within the medium. 

The remaining T = partial waves with L < 4 are both singlets: ^D2 and 
^64. However, the ^D2 phase shift, though positive over the energy domain 
of interest, is clearly dominated by the ^P2 phase shift, while the phase 
shift, again positive, is even less important. 

Therefore, in the conventional picture of neutron-star matter (see, however, 
ref. [63]), singlet and triplet pairing are operationally equivalent to ^Sq and 
^P2-^F2 pairing, respectively. There have been many illuminating studies of 
the former problem, notably [41,55,48,67,78]. Here we shall concentrate on 
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triplet pairing in the pure neutron system, ignoring the relatively minor mod- 
ifications that arise in beta-stable neutral matter with its dilute admixture of 
protons. 

BCS theory [105,106] furnishes a general framework for microscopic treat- 
ments of nucleonic pairing. As is the case for singlet pairing [55,67], two rather 
different problems must be overcome to realize a quantitative description of 
triplet pairing within this framework and make realistic predictions for the 
associated energy gap. The first problem to be solved is the construction of an 
accurate in- medium (or "effective") particle-particle interaction that provides 
the pairing force, together with a consistent mass operator that provides the 
single-particle input for the BCS equation. The second problem is to solve the 
BCS gap equation (actually a system of gap equations when one deals with 
J 7^ 0) , once the inputs for the pairing matrix and normal-state single-particle 
energies have been determined. 

This first task places great demands on even the most advanced tools of many- 
body theory. For ^Sq pairing, there exists a "first generation" of results on 
this facet of the problem, obtained by several different approaches (Landau 
Fermi-liquid theory [32], the method of correlated basis functions [41,55], the 
polarization-potential approach [47,56] and diagrammatic perturbation theory 
[66]). However, little progress has been made on this aspect of ^P2-'^F2 pairing, 
and the issue of modifications of the associated input pairing interaction and 
single-particle energies due to short- and long-range correlations inside the 
nuclear medium remains open. We shall not attempt to address this issue here. 
Rather, we shall be concerned with the task of actually solving the triplet BCS 
system, assuming the particle-particle interaction is identical with the bare or 
vacuum NN interaction, and that the input single-particle energies needed 
to complete the BCS equation(s) have the same form as for free particles 
(perhaps with an effective mass). This problem also proves to be nontrivial, 
since we face essentially the same difficulties as in singlet S-wave pairing, and 
more: 

(a) Again, the kernel of the BCS gap equation(s) has a pole at zero gap 
amplitude, implying an incipient logarithmic divergence. 

(b) Again, the nonlocality of the particle-particle interaction in momentum 
space implies slow convergence of the integral over momentum modulus 
k in the BCS equation(s) - precluding application of the standard weak- 
coupling estimates of BCS theory [106]. This feature of nucleonic pairing 
implies that one must integrate out to large values of the momentum 
variable. 

(c) As already intimated, an additional complication of pairing in higher 
angular momentum states with J 7^ arises from the (possible) angle- 
dependence of the gap, which can be parametrized in terms of depen- 
dences on the quantum numbers J, L of the total and relative orbital 
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angular momenta and on the magnetic quantum number M associated 
with J. One must deal with a set of coupled gap equations for the gap 
components A'l^[k) in an expansion of the (generally anisotropic) gap 



The combination of difficulties (a) and (b) can make numerical solution of the 
gap equation(s) a problematic exercise. To overcome these obstacles, a special 
procedure called the separation approach has been introduced in ref. [67] and 
applied to ^Sq pairing. In this approach, the generic pairing matrix V{k, k') is 
decomposed, identically, into a separable part, plus a remainder whose every 
element vanishes when either momentum variable lies on the Fermi surface. 
(As is characteristic of BCS theory, we assume that the particle-particle in- 
teraction is independent of frequency.) Applied to ^Sq pairing, this separation 
procedure reduces the BCS gap equation to two equivalent coupled equations: 
a nonsingular, quasilinear integral equation for the shape of the gap function 
and a nonlinear equation determining its amplitude at the Fermi surface. We 
hasten to point out that this strategy is by no means limited to the case 
where the pairing matrix elements are supplied by the bare NN interaction 
- it will retain its practical and analytical value when a realistic in-medium 
particle-particle interaction is employed. 

The primary aim of this article is to generalize the separation approach to 
pairing in an arbitrary angular momentum channel and apply it to '^P2 pairing 
in neutron matter. Since, according to consideration (c) above, we are now 
confronted - in general - with a set of coupled gap equations rather than 
a single one, the situation is much more complicated than for -^Sq pairing. 
However, in this case the power of the separation method becomes even more 
apparent. It enables a decomposition of the problem into (i) the determination 
of the shape of the gap function in /c-space through an appropriate nonsingular 
quasilinear integral equation whose kernel is practically independent of the 
energy gap, and, once this is done, (ii) the determination of a set of numerical 
coefficients specifying the angular dependence of the gap through a system of 
nonlinear 'algebraic' equations. 

In both singlet and triplet cases, the effectiveness of the separation scheme is 
greatly enhanced by the existence of a small parameter, given by the ratio of a 
suitable measure of the gap amplitude to the Fermi energy. Small-parameter 
expansions work well in the singlet problem, and even better in the triplet 
case. Corrections to the leading order in a small-parameter expansion of the 
sp^-Sp^ energy gap are down by a factor 10"^ or 10~^. This property leads to 
considerable simplifications in the analytical and numerical treatment of the 
triplet problem, greatly ameliorating complication (c). 

The separation approach is not merely a tool for numerical solution of gap 
equations. Following the lead of refs. [67,80,75,81], we shall demonstrate its 
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very important role in facilitating incisive analysis of the nature of pairing 
solutions in different angular momentum channels. In particular, the separa- 
tion formalism provides a basis for establishing certain universal properties of 
pairing, properties that are independent of density and independent of the 
specifics of the pairing interaction itself, whether furnished by the in- vacuum or 
in-medium particle-particle interaction. The analytical results obtained here 
and in refs. [75,81] remain valid outside the context of neutron matter and 
carry implications for the still more complex problem of pairing in liquid ^He 
[109]. 

In sect. 2, we display the BCS gap equations that describe pairing in an ar- 
bitrary two-body channel of given total spin S and isospin T. Sect. 3 extends 
the separation method of ref. [67] to this general problem, with special at- 
tention to ^P2-^F2 pairing. Sect. 4 is devoted to an analysis of pairing 
in the absence of tensor coupling to the ^F2 state ("pure" ^P2 pairing). We 
offer a straightforward criterion for determination of the upper critical density 
at which a solution of the pairing problem ceases to exist in the given two- 
body channel, in terms of vanishing of a characteristic determinant associated 
associated with the separation formulation. In sect. 5 we introduce the small- 
parameter expansion of the pairing problem and use it to clarify the nature 
of the angle-average approximation proposed by Baldo et al. [53]. In turn, the 
latter approximation is adopted to formulate an efficient procedure that al- 
lows reasonably accurate numerical solution of the coupled-channel, '^P2-'^F2 
pairing system when the goal is to find the magnitude of the pairing effect and 
the detailed structure of the pairing solutions is a lesser concern. We go on 
to recount how the fundamental character of the diverse solutions of the non- 
linear pairing problem may be elucidated within the separation method. The 
derivation of universal features of '^P2 pairing is exemplified and discussed. Re- 
sults from numerical application of the separation approach to triplet pairing 
in neutron matter are reported in sect. 6. Three representative models of the 
free-space NN interaction are used to construct pairing matrix elements, the 
Argonne Vis potential [107] being chosen as an example from the current gen- 
eration of models yielding high-precision fits of the A^A^ scattering data [108]. 
The results presented explicitly include the critical upper density for gap clo- 
sure, the momentum dependence of gap components, and the conventional 
average measure of the gap amplitude at the Fermi surface, variously calcu- 
lated for "pure" and tensor-coupled ^P2 pairing, with and without use of the 
angle-average approximation. These results are compared with previous find- 
ings of other authors. In sect. 7, we indicate some promising future directions 
for research on pairing in infinitely extended nucleonic systems. An appendix 
reviews the separation method as developed for ^Sq pairing. Small-parameter 
expansion is employed to devise an approximation for the gap amplitude that 
may be taken as a generalization of the BCS weak-coupling formula. Addi- 
tionally, the technique of small-parameter expansion is extended to the triplet 
pairing problem. 
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2 Any- channel gap equation 



To describe pairing in an arbitrary two-body state with total spin S, angular 
momentum J, and isospin T, it is useful to introduce a 2 x 2 gap matrix A(k) 
which has the expansion (cf. [57,30]) 

A(k) = Ai'\k)Gi^{n) (1) 

LJM 

in the spin-angle matrices 

Gi^(n)^ J2 {\lcT,a2\SMs){SLMsML\JM)YLMM: (2) 

MsMl 

with n = k/A;. The total spin and isospin quantum numbers S^T (= (1,1) 
for triplet states and neutron matter) are fixed throughout and generally sup- 



pressed. Setting (n) = Tr [G'^r*(n)G'^f'^ (n)J , the spin-angle matri- 

ces (2) may be shown to obey the orthogonality condition 

/ S'Jil'^^ (n)dn = ^vL^h-H^MM, . (3) 

The particle-particle interaction has the corresponding expansion 

y(k,k') = {k\VfL'\k')Gi''{n)Gi^*{n') , (4) 

LL'M 



with \L — L'\ < 2 in the case of tensor forces. The coupled set of generalized 
gap equations for the components A-[^{k) then reads 

aJiMi/t/n 



with A = S* + (L — L')/2 and dr' = (2/7r)(A;')^dA;'. The energy denominator 
-E(k) = [i'^{k) + D^(k)]^/^ involves the square of the single-particle excitation 
energy ^{k) = £{k) — /i of the normal system, measured relative the Fermi 
energy /x, and a term 

D\\.)^\ E ^i''*{k)Ai\^^{k)Sil'^''^{n) (6) 

^ LJMLiJiMi 



giving rise to an energy gap in the single-particle excitation spectrum of the 
pair-condensed system. Examination of the the explicit formula for the angular 
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factor S'/^'^^^^(n) reveals that a summation over total momentum Ji must 
be present in (5) and (6). (This fact is well known in the theory of superfluid 
^He [109]: the anisotropic ABM phase is built of states with total momentum 
J — 1 and 2.) The sum over orbital momenta L' on the l.h.s. of (5) appears 
only when tensor forces are present and is then restricted by |L — L'| < 2 and 
parity conservation. On the other hand, the sum over Li resulting from the 
angular dependence of the gap contains, in principle, an infinite number of 
terms. 

We identify the energy gap as |-D(k)| and note that it is in general angle de- 
pendent and in general has nodes. It is convenient to define an angle-averaged 
gap function A{k) through the positive root of 



A\k) ^ D\k) = i- / D\\.)dn =^\Y. \Ai^{k)\" . (7) 
AttJ 47r 2 ^jj^ I 



(The last equality follows from the orthogonality property (3).) Analogously to 
the ^So pairing problem, the quantity (7), evaluated at the Fermi momentum 
kp and denoted A^, = D'^ikp) furnishes an overall measure of the strength 
of pairing correction to the groimd-statc energy in the preferred state. The 
notation D'^ikp) is included here since it has been used earlier in the literature 
[57]. 

We take this opportunity to caution the reader that the literature is inconsis- 
tent and sometimes vague in the choice of energy gap measure. Care must be 
taken in interpreting statements about the magnitude of "the gap." Moreover, 
the choice made in normalizing or defining the pairing matrix elements is not 
uniform. We have attempted to be precise, and for the most part we follow 
the conventions estabhshed by Takatsuka and Tamagaki [24,57]. 



3 Separation formulation of the triplet pairing problem 



Generalizing the separation technique developed in ref. [67] for the ^Sq prob- 
lem, we again split the pairing matrix elements into a separable part and a 
remainder Wlj^,{k,k') that vanishes when either of the momentum variables 
/c, k' lies on the Fermi surface. Thus 

{k'\ViL'\k) = Vi^,{K k') = vi^^iuimiuik') + Wi^.iK k') , (8) 



where vi^ = V^L,{kF,kF) = V^/^^I/cf) and (pi^, = {kFlVLi'lkp) /vIl'- 
Substitution of this identity into the generic-channel gap equation (5) yields 
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= j:Di^..i>iAk) (9) 

where the coefficients 

Dif,-{-l)\L, UiAk) E SUt^^^^in) 2\ dndr(lO) 



(with dr = {2/7r)k'^dk) are just numbers. 

As in ref. [67] (see also the appendix), we have introduced a typical scale S 
in place of the gap in single-particle spectrum, in the term of (9) involving 
the portion Wli^,{k, k') of the pairing matrix that vanishes on the Fermi sur- 
face. Such terms are quite insensitive to the choice of this scale, since the 
vanishing of VF/^/(/c, k') at k' = kp ensures that the integral over k' receives 
its overwhelming contributions some distance from the Fermi surface where 
D'^{k') is negligible compared to i'^{k'). In the present context where the gap 
functions are angle-dependent, it is important to note that the 5 quantity may 
be regarded as angle-independent. Numerical calculations show that the final 
result is practically independent of the choice of 5 and we shall normally take 
it to be zero. 

Employing the orthogonality condition (3), one finds that only diagonal terms 
(meaning Li = L', Ji = J, M' — M) contribute to the sum on the l.h.s. of 
eq. (9), which becomes 

r A'^^(k'\ 
Ai^(A;)-i:(-l)^/ WUkX)^^dr' = Y.DiWLL^^^^ (11) 



Henceforth we shall frequently suppress the fixed index J. The functions 
4>LLi{k) = 4>^^^{k) superposed on the r.h.s. of (11) do not depend on M. 
Accordingly, the gap components A^(A;) = A'l^{k) can be constructed as 
linear combinations 

Af(A;)= E<l.Xl^^^(^) (12) 
L1L2 



of M— independent functions XiJ^^^ik) = Xl^^^^(^) multiplied by coefficients 
-^L^L2 determined by eq. (10). It is helpful to regard the quantities A|^(A;) 
as forming a vector with components indexed by L, and similarly the 
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functions XliLzI^) ^ composing a vector x^^^^. The equations 

Xi^'^ik) - El-l)^/ ^^^'(^'^O^l^dr' = S.L^cP'^^^ik) (13) 



for the shape factors Xi^^^ik) arc readily obtained from eq. (11). From this 
result we see that XL^^^ikp) = Slli for any L2, since the potentials WLL'{k, k') 
vanish for k on the Fermi surface while ^^^^^(A;^) = 1 for any Li, L2. It is im- 
portant to observe that the equations (13) are free from any log singularities 
and form a coupled system of linear integral equations. We further remark 
that these equations can be solved without regard to the values of the coeffi- 
cients Di^, = D^L, given by the set of equations (10) with Ai\^'{k) = Af^ 
expanded as in (12). 

It is also important to take note of an essential difference in the reformulated 
equations (11)-(13) as compared to the original gap equations (5). The sum 
over orbital angular momenta Li in (5) is in principle infinite, while in (11) 
only the states with V — L and L' — L±2 can contribute, due to the conser- 
vation of the total angular momentum J and parity. Therefore the number of 
components of the vectors A^^ and x^'^^'^ docs not exceed 3, since \L~ L'\ < 2 
and in eqs. (11)-(13) the total angular momentum J has a fixed value. We 
may infer that, in the case of a ^P2-^F2 mixture where J — 2, orbital angu- 
lar momentum quantum numbers beyond 5 contribute only to the numerical 
factors Dff^^^. 

It is instructive to write out eqs. (11) and (13) in extended form for the 
^P2-^F2 problem. Using a symbolic notation in terms of the operator T — 
Wl2\i^ + (52]V^ the system (11) reads 

Af + Tn Af + Ti3 Af = L>^(0i^ + r>j^0^=^ , 

A3^ + T31 Af + T33 Af = L>3^0=^^ + L>3^(/>^^ . (14) 

Suppose all the coefficients D^^^i are zero except D\[ ; then the solution A^ 
of (14), a two-component vector, reduces to 

Af(^) = <xi^(^) (^ = 1,3), (15) 
where the components x^ and x^ of xx\ obey the pair of equations 

X3'+r3lXi'+T33X3' = 0. (16) 

The vector xis is given by a second pair of equations. 
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xf + TsiX? + Tssxf-0, (17) 

derived by supposing that is present but all other Df^'j^, are zero. The 
remaining two pairs of equations for the vector functions and x^^, namely 



Xf+Tsixf+Ts^xf 



k31 



(18) 



and 



xf + Tnxf + Tuxf^O, 

xf + T3,xf + n,xf (19) 

are constructed in the same fashion. The general solution of the ^P2-^F2 pair- 
ing problem is given by eq. (12) with the indices Li and L2 taking on the 
values 1 or 3. It is worth noting that the same system of equations arises 
when one deals with S — D pairing in the deuteron channel, except that the 
indices Li and L2 then take on the values and 2. 

Let us return now to the problem of pairing in an arbitrary two-body channel. 
Having learned that the shape factors XL^^^{k) in the representation (12) may 
be found by solving a set of nonsingular linear integral equations independently 
of the coefficients Dff^^^, we must next consider the determination of these 
coefficients, assuming knowledge of the x functions. Substitution of (12) into 
(10) provides a set of nonlinear 'algebraic' equations to be solved for the 
desired set of numbers: 



)JiMi 
L2L3 

Li J1M1L2-L3 

JL2LS I 



X /^i' '(^^ dndr. (20) 



In the vicinity of the (upper) critical value Pc of the density p where the gap 
value vanishes, the term -D^(k) in the denominator of (20) can be omitted, 
and in this case we arrive at a system of linear equations for the coefficients 
D'l^,. The determinant of the linearized system should vanish at p = Pc, 
since by definition there is no solution at that point. Thus we obtain a new 
condition for gap closure. For example, in the case of ^P2-''F2 pairing (where 
channel coupling is included), the zero of the determinant corresponding to 
the dominant V{{^'^ matrix elements of the Argonne Vis potential [107] gives a 
value kc — S-Q fm~^ for the critical Fermi momentum when the input function 
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^{k) — e{k) — // in eqs. (20) is constructed from free single-particle energies 
e{k). 



In the standard numerical approach to BCS pairing problems, the original 
system (5) for the gap components A£^(A;) is solved by iteration. One starts 
with suitable "guesses" for these functions (ten in number, for ^P2-^F2 pair- 
ing) and inserts them into the r.h.s. of (5), generating a new set of component 
functions; and so on until convergence is hopefully achieved. Such a procedure 
encounters obstacles even in the case of ^Sq pairing [67,48,55]. A conventional 
iterative approach is even more problematic for pairing in higher angular mo- 
mentum states, where the dimensionality of the nonlinear problem is much 
larger. From this point of view, Eqs. (20) are far more appropriate as basis 
for numerical computation, since shape factors XL^^^i^) are determined sep- 
arately and iterations are needed only to find the set of coefficients -D^*^^ 
entering the representation (12). 

In principle, instead of trying to solve the system (20) directly, we could use 
the variables Dj^^, (twenty in number, for ^P2-^F2 pairing) to parametrize 
the condensation energy and apply the variational principle to determine an 
appropriate set of these factors and a corresponding set of gap components 
^jj^ {k). However, one must remember that the performance of optimization 
algorithms, as well as the convergence of iterative procedures, is very sensi- 
tive to the degeneracies of the functional under consideration. Accordingly, 
any attempt at implementing such an approach should be prefaced by a thor- 
ough study of the symmetries of the gap problem and the degeneracies of its 
solutions. Such a study is of considerable fundamental interest in any case 



4 Triplet-P pairing without channel coupling 

To gain a better understanding of the triplet pairing problem in neutron mat- 
ter, we revisit the pure (i.e., uncoupled) ^P2 case analyzed in some detail in 



If tensor coupling to the ^F2 state is neglected, the system of gap equations (5) 
reduces to a set of five equations, corresponding to M = 0, ±1, ±2. Taking into 
account the property of time-reversal invariance, which implies the relations 



it is seen that only three of the generally complex components are in- 
dependent. The independent components have customarily been written as 



[75]. 



[24,57,40]. 




(21) 
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[24,57,40] 



Af{k) = {Al'-\k)y = 62{k) + t7j2{k) , (22) 
where the Si{k) and rjilk) are real functions of k. 

Takatsuka and Tamagaki [24,57], as well as Amundsen and 0stgaard [40], have 
carried out numerical studies of the pure problem that demonstrate an 
interesting feature of triplet pairing in neutron matter (and in fact of triplet-P 
pairing more generally). Solutions of different types were found to be nearly 
degenerate, i.e., they yield very nearly the same condensation energy. The 
spread of condensation energies reported by Takatsuka and Tamagaki at a 
typical density is only 3%, and 0stgaard and Amundsen obtained very similar 
results. The calculations of these authors were based on an older nucleon- 
nucleon potential called OPEG (see sec. 6 for details). Although this potential 
model is no longer competitive, its essential aspects are semi- quantitatively 
correct. At any rate, it will be argued that the nearly degenerate character of 
the level structure in the triplet pairing problem and the energy ordering of 
different solution types are very insensitive to the interaction assumed. 

Important insights into the uncoupled problem can be achieved within the 
separation method [75]. For pure ^P2 pairing, eq. (15) becomes simply 

Al''{k)^Dl^xl\k), (23) 
while the shape factor xi^ (^) obeys the uncoupled equation 



We now make the decomposition D^2^ = (Am + ifiM)So{kF) / where 
Soikp) — Df(2^ serves as a scale factor. Relations between the parameters 

Ajv/ and Km and the quantities Si{kp) and rii{kp) {i — 1,2) stem from the 
definitions contained in eqs. (22). Evaluating the spin sums 

2M' and sep- 
arating real and imaginary parts in eq. (20), we extract the following set of 
equations determining the coefficients D{f'^ (cf. refs. [75,81]): 



A2 = -v'[\2{Jq + J5) - Ai Ji - K1J2 - J3] , 
«2 = -v'\k2{Jq + J5) - i^\J\ + Al J2 + J4] , 
Ai = -i;'[Ai Je - (A2 + 1) Ji + ^2 J2 - /«i J4/2] , 
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= —v'IkiJ'j — K2J1 — (A2 - 1) J2 - Ai J4/2] , 

1 = _^;'[_(Ai Ji - J2 + A2 J3 - «2 J4)/3 + J5] , (25) 

where v' = (7r/2)f/j"^. The integrals Ji are given by [75] 

, fr^n , 0^^(fc)xP(fe) .2..dn 



with /o = 1 - 3^2^ /i = 3x^/2, /2 = 3y^/2, /3 = 3(2^2 + - l)/2, /4 = 3xy, 
and /s = (1 + 32;^)/2, where x = sin cos y = sin^^sin(/9, and 2; = cos^^. Of 
these integrals, only J5 contains a principal term behaving like InAp, since 
only /5 among the fi does not vanish under the angular integration when 
multiplied by an angle-independent factor. 

The normal-state single-particle energy ^(k) present in the energy denomina- 
tor E(k) = [^^(/c) + Z^^(k)]^/^ is conventionally assumed to have finite slope 
d^{k)/dk, while the gap term takes the explicit form 



D'{k) = -^[(1 + A2)^ + + 4 + (A? - 4A2 - kI)x' 

- 2(Ai + A1A2 + KiK2)xz + (3 + A^ - A^ - 2X2)z'^ 

+ 2(2^2 - KiXi)xy + 2{ki + X1K2 - X2Ki)yz] [xl\k)]^ ■ (27) 



In their angular content, the system (25) and the trigonometric decomposi- 
tion (27) are consistent with the formulations of the uncoupled problem 
given in rcfs. [57,24,40], and specifically with eqs. (4.7)-(4.8) of Amundsen and 
0stgaard. However, our separation formulation, as expressed in eqs. (23)~(27), 
has distinct advantages as a basis for numerical solution of this problem. To 
begin with, there is only one shape function x{k) = Xp(^), which satisfies a 
linear integral equation of exactly the same form as for singlet-S pairing (cf. 
eq. (60). This property of shape invariance follows naturally from the formal 
development presented in sect. 3. Given the key physical fact that the an- 
gle dependence of the problem only comes into play significantly in the close 
vicinity of the Fermi surface, we observe that this part of momentum space is 
almost entirely eliminated from the k' integration in eq. (24) by the vanishing 
of W{k, k') at k' = kp. Thus, to an excellent approximation, all solutions Af^ 
of the uncoupled triplet problem have identical shapes. In particular, the po- 
sitions of their zeros and maxima coincide, as is clearly seen in the figures of 
previous works on triplet pairing in neutron matter [57,40]. 

To establish the angular dependences of the allowed solutions of the •^P2 prob- 
lem, the shape equation (24) must of course be supplemented by a set of 
equations that determine the coefficients or amplitudes corresponding to the 
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different M values. An important benefit of tfie separation strategy is tliat 
the nonlinear aspects of the problem are now concentrated entirely in the lat- 
ter system of equations, given by (25) - equations for a set of numbers 
rather than for a set of functions. By contrast, earlier authors [24,57,40] have 
worked directly with a set of singular linear integral equations for the five 
functions Soik), 5i{k), rii{k), 52{k), and r]2{k) defined in eq. (22). In this ap- 
proach, the dependences of gap quantities on the modulus of the momentum 
and on angles (or magnetic quantum numbers) are confronted simultaneously, 
with attendant opportunities for numerical difficulties. 

Obviously, the process of numerical solution is facilitated by the separation 
of "shape" and "angular" aspects of the problem that has been accomplished 
in eqs. (23)-(27). Less obviously, this transformation opens the way to an 
incisive analysis of the nature of these solutions. Such an analysis has been 
carried out in refs. [75,81], extending the formalism to finite temperature T 
via introduction of the factor tanh [£'(k)/A:BT] in the integrands of the basic 
gap equations eq. (5) and the Ji integrals (26). 

Appeahng to the orthogonahty of the spherical functions for different M, it is 
possible to find special solutions of the system (25) (or of the system (4.7) of 
ref. [40]) that satisfy most or some of its equations identically. The simplest 
types of consistent solutions are characterized as follows [24,57,40]: 

1. Only M = ±2 components are present: So = Si = rji — 0. 

2. Only M = is present: Si = rji = S2 = f]2 = 0. 

3. Only M = ±1 are present: Sq = S2 = TI2 = 0. 

4. Coupled case containing M = 0, ±2: Si = rji = 0. 

In addition, there is naturally the general case involving all five M values. 

For purposes of illustration, let us examine the situation for solution Types 1 
and 2 in some detail. 

Type 1. Under the assumption that only the component with M = 2 is present 
(its M = —2 partner being included trivially), we need only deal with Ap(A;) = 
S2{kF)x{k): noting that a free choice of phase allows us to set 772 = 0. In this 
case only the first of eqs. (25) is relevant and it becomes 

1 f , . , s i sin^ 9 , , 

- = - / (f)(k)x{k)^^ ^ dndr , 28 

^ -f 2^e'(A;) + |A|x'(A;)sin2^ 

where v = v(j"'^. All the other equations of the set (25) are met identically. 
(We remind the reader that is the value, on the Fermi surface, of the 
angle average defined in eq. (7).) Apart from the angular dependences and 
angular integration, eq. (28) has the same form as the equation (61) for the 
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gap amplitude in the ^Sq problem (see the appendix). The angular integration 
can actually be performed analytically, with the result 



- = -/ 

V J 



mx{k) 



$2(/32)dT, 



(29) 



where 



Aarcsin/3. + A 



1-0: 



2 arcsin P2 



02 



(30) 



The root of eq. (29) can be found numerically without difficulty, after find- 
ing the universal shape function x{k) = xV{k) by matrix inversion of the 
discretized form of eq. (24). This has been done for pairing matrix elements 
calculated from the OPEG potential employed by Takatsuka and Tamagaki 
[24,57] and Amundsen and 0stgaard [40]; some results for this case are re- 
ported in sect. 6. 

Type 2. Here there is the single component Af{k) = So{kF)x{k) = Apxik)- 
The dispersion equation determining Apis given by the last of eqs. (25), which 
becomes 



jmx{k)- 

V J ' 



i(l + 3cos2e)) 



2Je{k) + '.Al{l + 3cos^e)x'{k) 



=dndr 



(31) 



Analytic integration over angles is again possible, yielding 

mx{k) 



1 

V 



2^e(k) + lAlx'ik) 



M0o)dT, 



(32) 



where 



(33) 



iPoik)) 



e{k) + \Aixm 



Some numerical results for this case are also reported in sect. 6, again for the 
OPEG potential. 
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5 Small-pcirameter expansion and angle-average approximation 

Quantitative description of pairing in nuclear systems is beset with special dif- 
ficulties because the strong inner repulsion (or strong momentum dependence) 

present in some partial waves (notably S states) implies that the momen- 
tum integral entering the gap equation must be extended to large momenta 
[48,55,67]. This behavior precludes the application of simple estimates such as 
the BCS weak-coupling formula, which assumes sharp localization of the pair- 
ing interaction around the Fermi surface. However, important simplifications 
are permitted by the existence of a small dimensionless parameter dp-, given 
essentially by the scale of the energy gap at the Fermi surface, relative to the 
value of the Fermi energy Sp- 

In the case of ^Sq pairing, the small parameter is just dp = ^f/^f-, where 
Ai? = A{kF). Framing this problem within the separation method, Taylor- 
series expansion of the energy denominator of the equation for the gap ampli- 
tude Ap has been employed in ref. [80] to derive a sequence of approximants for 
this quantity expressed as series of terms in (dp)'" In dp and (dp)'^ with integral 
n. The derivations involved are outlined in appendix A. Even the leading (ze- 
roth) approximant yields reasonable results: the exact solution for Ap within 
the relevant density range is reproduced within 5%, based on pairing matrix 
elements formed from the ^Sq component of the Argonne Vis interaction. 

In the triplet case considered here, the small parameter is appropriately defined 
as dp = [D'^(kp)Y/'^ /ep = Ap/sp. A characteristic value of this quantity in 
the density region of interest is 0.5 MeV/80 MeV ~ 6 x 10~^, which is at least 
an order of magnitude smaller than typical values for singlet-S pairing. On this 
basis, one may expect to achieve good accuracy through an approximation that 
retains only the principal logarithmic terms ~ InZ)^(k) together with those 
independent of D^(k), while omitting corrections ~ d'jj^lndp. Small-parameter 
approximations are examined in more detail in the appendix. Here we restrict 
ourselves only to brief remarks. 

Consider the integral on the r.h.s. of eqs. (20) for the coefficients Df^. We 
may carry out an integration by parts with respect to the variable ^, to derive 
a principal term (within the angular integration) that goes like lnD'^{kp,n) 
and a remainder that involves -D^(k) only in the factor 

ln{\C\+[e + D\kf'). 

Since latter quantity is quite smooth in the vicinity of the Fermi surface, its 
replacement by ln(2|,^|) is valid within corrections of order of d^plndp that 
belong to the next stage of approximation. Writing D^{kp,n) = ApV^{n), 
we extract the scale Ap and express the angular dependence by I>^(n). The 
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dominant log term then decomposes as lnD^(/ci?,n) = InA^, + lnl?^(n). By 
virtue of the orthogonahty condition (3), the angle-independent portion In A^^ 
contributes only to terms on the r.h.s. of eq. (20) that are diagonal in the 
magnetic quantum number, i.e., only for Mi ~ M. Since is a small quan- 
tity, the contribution from the scale term In A% greatly exceeds that from the 
angle-dependent part Inl)^(n), which remains of order of 1 even in the limit 
A^ ^ 0. It is worth noting that the A-independent contributions coming 
from terms in ln(2|^|) are also of order of 1. However, these terms are angle- 
independent and, like the scale term, contribute only to the diagonal part of 
the r.h.s. of (20). In fact, to within corrections of order ~ dplndp, the result 
for the D coefficients calculated from the angle- independent terms alone (i.e., 
dropping InD^(n)), is the same as the result we obtain if we simply replace 
the quantity -D^(k) in the energy denominator of the original equation by its 
angle average (7), 

^(^) = ^EiAr^°(^)r +^ E iAn^)r- (34) 

l LM>0 



We reiterate that justification of this replacement, which corresponds to the 
angle-average approximation iiist introduced by Baldo et al. [53], rests on (i) 
the unimportance of the angle-dependent quantity InP^(n) relative to In A^^ 
and (ii) the unimportance of correction terms of order d^rlridp or higher in 
the small parameter expansion. The angle-average approximation generally 
overestimates the pairing gap and the magnitude of the pairing energy (i.e. 
the pair condensation energy correction to the energy of the normal ground 
state). However, it allows a considerable simplification of the triplet pairing 
problem with very modest sacrifice in accuracy, if our primary interest lies 
in the magnitude of the pairing effect rather than the detailed structure of 
pairing solutions. 

Making the replacement (34) in the general coupled-channel equations (20) 
for the gap amplitudes, we arrive at the quasi-uniform system 

D¥l' - (-1)^..' / E <L.-r^^=^^dr . (35) 



Since the index M appears explicitly only on the quantities being sought, there 
exists a degeneracy with respect to this quantum number: All sets of -D^ff 
yielding the same result for the gap quantity D'^{k) of eq. (34) are equivalent. 

We first apply the angle- averaging prescription to the pure problem. In 
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this case eq. (35) becomes simply 



l = -.I^M^ir (36) 



with D'^ik) = Ap[x\'^{k)]'^ = Aj^x^ik) and v = vf^"^. This equation, having 
no explicit reference to magnetic substatcs, is identical in form to the singlet- 
channel equation (61) and is solved by the same method. Numerical results 
from solution of the uncoupled gap problem in this approximation are pre- 
sented and discussed in sect. 6. If we do not care about the detailed angular 
behavior of solutions, but only about the size of the gap (as measured by Ap) 
and the condensation energy, the angular-averaging prescription gives very 
good predictions. This statement is justified by the close agreement with the 
results shown in fig. 1 for solution Types 1 and 2 and by the near degeneracy 
of the Ap values and condensation energies found numerically [24,57,40] and 
analytically [75] for the different solution types. 

Turning to the general coupled-channel case, observe that the terms with 
Li ^ V on the r.h.s. of (35) contain integrals 

0^^'(A;)xi/^^(A;) / 0^^'(^)xi/'^^(A:) 



2Je{k) + D\k) J 2|e(A;)| 



J mk)\ — ^^'^ 



which are almost independent of the gap value, since the property Xi)^^ — ^lli 
imphes that the numerators of the integrands vanish at the Fermi surface. In 
the diagonal case, Li — L', we introduce an unknown positive number 

-Ar (38) 



and write 



0^^'(fc)x^:^^(fc) 



^T = Z^Blt,\ (39) 



where 



^ / 0^^'(fe)x^:^-(fc)-[xi^w]^ ^^ ^ s 0^^'(fe)xi:^-(fe)-[xiHfe)]^ ^^ 

^ 2Je{k) + D^{k) ^' 2|e(A;)| 



The simplification on the right in eq. (40) is permitted since (j)^^'{kp) ~ 
Xu^^ikp) = 1 and the numerator in the integrands again vanishes at the 
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Fermi surface. The integrals -Bf^^ are effectively already known: with the 
indicated simplifications, they only depend on the x functions, which may 
be found from eqs. (13) without reference to the gap quantities D^y being 
sought. 

We now specialize to coupled-channel, ^P2-^F2 pairing where the system (35) 
for the coefficients Df^j^, contains four equations. For a solution of this set of 
equations to exist, the determinant 



C{Z)^ 



(l/vn) + Bll + Z 

-(iMs) 



^13 

Bll + Z 



/1 13 
^13 



^33 



BM + z 

431 

^33 



All 

Bfl + Z 

-(iMi) + Ail 



^33 



7^33 



^31 



,(41) 



Bll + Z (1/^33) + Bll + Z 



which depends on Aj? (or the -D^j,/ quantities) only through the parameter Z, 
must vanish. Evaluation of this determinant may be simplified by subtracting 
the third row from the second row and the second row from the fourth row, 
and then subtracting the second column from the first column and the third 
column from the fourth column. This leads to a quadratic equation for the 
determination of the unknown parameter Z: 



C{Z) = C2Z^ + ciZ + Co = , 



(42) 



where 



C2 



Co = C(0) 



C'oo C'oi ^"02 Co3 

Cio Cu C12 Ci3 

C20 C21 C22 C23 

C30 C31 C32 C33 



C30 ~ CiQ — [C31 — Cii] C33 — Ci3 — [C32 — C12] 

Coo ~ C*20 — [C'oi — C21] Co3 — C23 — [C02 — C22] 
Ci = C(l) - Co - C2 . 



(43) 



(44) 
(45) 



Equation (42) can have two roots {Zq}, and only positive solutions are of 
interest. 

Having found the appropriate root of (42), one may solve any three of the 
four equations of the system (35) for the ratios Dffj^,/D^f. These ratios may 
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then be used to obtain an equation for that is analogous to the algebraic 
equation for the gap amplitude Ap in the ^Sq problem, namely 

I^"<*)1" -Ar. (46) 



2Je(.k) + D'(k) 



Solution of this equation yields the desired gap function in angle-average ap- 
proximation. 

The angle-average approximation serves well if the goal is to calculate the gap 
amplitude Ap = [/^^(fci?)]^/^. However, if information is required on the spec- 
trum of pairing solutions and the energy splittings between these solutions, it 
is essential to include the angle-dependent effects from the terms in eq. (20) 
going like InP^(n). By so doing, accurate results and valuable insights can be 
gained even if we stop at leading order in the small-parameter approximation. 

The situation is made more tangible by considering two simple solutions in- 
volving only one M value (and its negative), specifically M = 2 or M = 
(corresponding to Types 1 and 2, respectively). Working at lowest order in 
the small-parameter expansion, partial integration of the definition (26) of 
the integrals over the variable ^ leads to the simplified form 

Ji--2l I Mz,^)lnD'{z,ip)^, (47) 





having made use of the equalities (f){^ = 0) = 
eqs. (28) and (31) are recast, respectively, as 



0) = 1. With this result. 



1 = -^J5 (Af =^)) and 



1 = 



Jo (Ar°)) + J, (Ar°))] . (48) 



The splitting between the states in question is determined by subtracting one 
equation from the other. Straightforward manipulations give 

/a(')\' 2 2n 



When the angle-average approximation is applied, eq. (28) or (31) is replaced 
by 



+ / de + / de — + In -2- - In — = , (50) 



vk^p J 2e J 2|e| A| 87r 



-1 
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which is to be solved for Ap = A'^. (Here, (p = 0(e) and x = after a 
change of integration variable from A; to e = {k/kp)^ — l.) Evaluating the ratios 
of A^"* and A^'* to A^^ along the lines of the above derivation of A^p /Ap\ 
the ordering of the three gap amplitudes is seen to be A^-* < A^'' < A^. 
This order agrees with the numerical results appearing in fig. 1 and with the 
findings of other authors [24,57,40]. 

It is interesting and highly significant that the result (49) for the ratio of pair- 
ing energies for the Type 1 (|M| = 0) and Type 2 (M = 0) solutions - obtained 
to leading order in the small parameter of the theory yet very accurate — is 
completely independent of the parameters of the actual interaction and inde- 
pendent of density. In other words, it is universal. 

Proceeding to the next layer of complexity, let us consider the Type 4 solution, 
again studied in leading order. In this case, M — and \M\ — 2 substates 
coexist, and the system of equations that determines the coefficients Df^ 
becomes 



A2 = -V [A2(Jo + J5) - -^3] , 
1 = -V [-A2 J3/3 + J5] , 



(51) 
(52) 



where we have introduced the more compact notation A2 = S2VG/S0 and the 
Ji are given by (26). The first of these equations can be conveniently rewritten 
as 

= A2Jo + (A^/3-l)J3 (53) 



by subtracting A2 times eq. (52) from eq. (51). 

The alternative form (53) is distinctly advantageous: unlike J5, the integrals 
Jo and J3 do not contain a term proportional to lnA|. and hence are inde- 
pendent of any details of interaction. Consequently, eq. (53) provides a closed 
equation for the parameter A2 whose solution is universal - i.e., the same for 
all pairing interactions. Once the universal parameter A2 is determined, the 
scale parameter 60 is calculated from (52); its value does, of course, depend 
on the interaction. 

The formulas for Jq and J3 can be further simplified by partial integration 
of Jo over z and J3 over (f. The remaining integration over (p is performed 
analytically. After simple operations, we arrive at 

1 \ r^H^ 3 -I- 3 }■ 

J„(0,A.)=--+4/ J3(0.A.) = ^- //Wd.. (54) 
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where f(z) = [(1 + + (3 - Xl)z^y - 24Ai(l - ^2)2)^2 general, these in- 
tegrals are expressible in elliptic functions, but in exceptional cases they reduce 
to elementary integrals. This is exactly what happens when we seek roots of the 
system (51-(53). In addition to the root A2 = corresponding to the M = 
state (the Type 2 solution), we find the roots A2 = ±3 and A2 = ±1. It should 
not be surprising that the roots are degenerate with respect to signs, since 
eq. (53) is invariant under reversal of the sign of A2. What is most surprising 
is that these roots are rational numbers. 

Rationality of the roots A2 is a consequence of a hidden symmetry of the prob- 
lem that transcends the approximation we have employed. To substantiate this 

assertion, consider first the particular solution A2 = 3, for which the squared 
gap function (27) becomes D'^{k, x, z; \2 = 3) ~ [xi^(/i;)]2[4 - 3(a;2 + z"^)]. 
The symmetry of this function with respect to x and z implies the rela- 
tion 3Jo(0, 3) -I- 273(0, 3) = 0, since the combination 3/o -|- 2/3 = 6{x'^ — z'^) 
changes sign on interchange of x and z whereas and other factors within 
the integrand of (26) are left unchanged. The above relation between Jq and 
J3 satisfies eq. (53). The essential point is that this relation derives from a 
simple symmetry of the problem; it is not predicated on any approximate 
scheme for the evaluation of the Jj. Now consider the root Ai = 0, for which 
D'^(k,x,y, z) ~ [xi''^(^)]^(l — x"^)- Since D'^{k,x,y, z) is independent of y and 
z, one may trivially integrate (26) over these variables. Simple algebra then 
gives 3Jo(0, 1) = 273(0, 1), and this relation also satisfies Eq. (53). We may 
therefore conclude that the ratio of M = 2 and M = components in the 
Type 4 solution is independent of any input parameters including density, 
temperature, and the parameters specifying the pairing interaction. 

A general analysis along similar lines, again exploiting the simplifications pro- 
duced by the separation method, reveals the fundamental structure of the full 
set of solutions of the BCS ^P2 pairing problem in neutron matter [75] . Con- 
sidering general multicomponent solutions, the allowed ratios of coefficients of 
different spin- angle matrices G"[^^{n) in the expansion (1) of the gap matrix 
turn out to be universal. As in the above derivation, any reference to inter- 
action, density, or temperature cancels out to very high accuracy in forming 
the ratios. The complete set of solutions for these ratios is obtained exphcitly 
by solving a set of five coupled 'algebraic' equations, one of which serves to 
determine the energetic scale of the pairing effect. The solution set and the 
spectrum of pairing energies are found to be highly degenerate in character. 
The small splittings between different solutions can be found by manipulations 
similar to those performed above for the Type 1 and Type 2 solutions. 

In principle, analogous analytic arguments may be used to establish universal- 
ities of pairing in the much more difficult problem of liquid '^He, where states 
with L = S = 1 and J = 0, 1, 2 contribute on an equal footing and the number 
of equations to be solved rises from five to nine. 
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6 Numerical applications: uncoupled and coupled channels 

The separation approach to BCS-type pairing problems which has been further 
developed and studied in this paper leads to considerable simplifications of the 
process of solving the gap equations in a generic pairing channel. Aided by the 
existence of a small parameter, this approach allows one to find the interval of 
existence of solutions without having to perform any iterations, and it allows 
one actually to find solutions without having to iterate a set of integral equa- 
tions. Some iterations may be required to construct solutions, but they need be 
performed only on a set of algebraic equations for the amplitudes of the com- 
ponents of the gap matrix (quantities that are momentum-independent but 
dependent on the magnetic projection M). The "division of labor" that un- 
derlies the separation method thus facilitates the process of finding solutions; 
moreover, because of the reduced number of numerical operations, it permits 
higher accuracy to be achieved. In addition to these numerical advantages, the 
new method offers an efficient framework for fundamental analytic investiga- 
tions of pairing, as has been demonstrated in sect. 4 and more systematically 
in refs. [67,75,80,81]. 

A necessary step in testing our procedures for numerical solution of the triplet 
pairing problem in neutron matter is comparison with results obtained by 
other authors. As indicated in sect. 4, Takatsuka and Tamagaki [24] (see also 
their review [57]) have studied the pure pairing problem at some depth, 
classifying the types of solutions and generating numerical solutions by iter- 
ation. Subsequently, Takatsuka [30] (sec also [57]) implemented an iterative 
procedure for solution of the gap equations in the '^P2-'^F2 coupled-channel 
case with fixed M = 0. This early numerical work was based on several older 
potential models [110] with soft (Gaussian) cores. The Takatsuka- Tamagaki 
calculations were repeated and supplemented by Amundsen and 0stgaard [40] . 
Baldo et al. [53] introduced the idea of angular averaging in considering M — 
and \M\ = 2 solutions (i.e., Type 2 and Type 1 solutions) of the pure case 
and reported preliminary numerical results for a separable approximation to 
the Argonne Vi^ potential [111]. More recently, Elgar0y et al. [69] have per- 
formed uncoupled and coupled ^P2-^F2 calculations with fixed \M\ (0 or 
2) in /3— stable matter for the Bonn potential models designated A, B, and C 
[112,113]. 

In the most recent microscopic study of triplet pairing, Baldo et al. [79] ob- 
tained results for several realistic models of the nucleon-nucleon interaction, 
including four that may legitimately be regarded as phase- shift equivalent. 
The latter potentials - Argonne Vis [107], CD-Bonn [114], and Nijmegen 1 and 
11 [115] - belong to a new generation of models that yield high-precision fits 
of the sanitized world supply of pp and up scattering data up to a bombard- 
ing energy of 350 MeV. All four produce fits with per data point close to 
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unity. In addition, gap calculations were performed for three "older" poten- 
tials - Argonne U14, Paris [116], and Bonn B. The authors of this extensive 
study conclude that reliable predictions of the ^P2-^F2 pairing gap at densities 
above about 1.7 times the saturation density of symmetrical nuclear matter 
will require potential models that fit the nucleon-nucleon phase shifts up to 1 
GeV. 

For the purpose of comparison with the pre-1990 results on the triplet pairing 
problem, we have chosen the OPEG ( "one-pion-exchange Gaussian") poten- 
tial [110] used in refs. [24,30,57,40], and specifically the version denoted OPEG 
^0 — 1. This potential has been fitted to scattering data collected before 1970, 
and it cannot be regarded as a realistic model of the nucleon-nucleon inter- 
action. However, it possesses all of the essential qualitative features of the 
interaction and describes them at a semi-quantitative level. Of local char- 
acter, it contains central, tensor, spin-orbit (L • S), and quadratic spin-orbit 
(e.g. (L- S)^) components, as written out in detail by Amundsen and Ostgaard 
[40]. 

To compare with the latest calculations [79] , we have generated results for the 
Argonne vu and Argonne Vig models. The diagonal pairing matrix elements 
of these two potentials are plotted in fig. 2. 

In presenting the numerical results, we employ two measures of the gap ma- 

r -,1/2 

trix, namely (i) the average quantity Ap = D'^{kp) defined via eq. (7) 
or eq. (34) and (ii) the gap component AjJ^{k), either as a function of the 
momentum modulus k or evaluated at the Fermi momentum kp- We note that 
when the pairing problem is solved in the angle-average approximation speci- 
fied in sect. 5, the M dependence of the latter measure becomes moot. Thus 
the index M will often be suppressed, along with the fixed index J = 2. 

The zeros of the characteristic determinant of the linearized version of eq. (20) 
provide us with candidate values for the upper critical Fermi momentum k^. 
This is true whether we restrict ourselves to pure ^P2 pairing or address the 
coupled-channel '^P2-'^F2 case. Examining first the uncoupled case, in fig. 3 we 
superimpose plots of the diagonal matrix elements = V/f"'^{kF, kp) of the 
Argonne Vi^ potential, the associated characteristic determinant, and the gap 
component Al^^kp) = Ai^kp) obtained in the angle-average approximation, 
all considered as functions of the Fermi momentum kp- From the preceding 
formal developments it is evident that a zero of the potential implies a pole of 
(j)^^{k) and hence a simple pole of the determinant. Such a pole is very promi- 
nent in fig. 3 and is naturally followed closely by a zero of the determinant. 
The position of this zero gives a value k^ — 3.6 fm^^ for the upper critical kp 
value at which the gap closes, the corresponding upper critical density being 
Pc = Oil the plot, the gap amplitude is clearly negligible at this point. 

Because of the exponential falloff, it appears to be zero even below this point. 
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but is in fact still finite over an interval in kp where vd ^ has positive values. 

Fig. 4 displays the characteristic determinant in the coupled-channel case, 
again with the Argonne f is potential as input for the pairing force. The deter- 
minant is evaluated for three different effective-mass ratios m* — M*/Mm the 
free normal-state single-particle spectrum €{k) = fi^k"^ /2Mn being replaced by 
h /c^/2M* to simulate dispersive effects of the medium. For any of the chosen 
values of m*, the determinant has a pole near kp = 1.9 fm~^ followed closely 
by a zero near kp — 2.2 fm~^, and a second pole near kp — 3.3 fm~^ followed 
closely by a zero near kp — 3.5 fm~^. It is the second zero that should be 
identified with the upper critical density, since the adjacent pole corresponds 
to a zero of f j'f ^. The point is that we are seeking to describe T = 1 triplet-P- 
wave pairing as modified by the tensor coupling to the T = 1 triplet-F channel. 
The first pole of the determinant corresponds to the zero of ^33""^- Thus the 
accompanying determinantal zero is not relevant to the task at hand, which is 
the identification of the upper critical point for triplet-P rather than triplet-F 
pairing. Generally, a lower hound on k^ for pairing predominantly in a state 
with orbital angular momentum quantum number L is given by the location of 
the pole that is produced by the vanishing of vlj^ (we assume that this occurs 
only at one density). Conventional wisdom gives < as a sufficient (but 
not necessary!) condition for a pairing instability in partial wave L. The in- 
dicated lower bound property is of course consistent with this statement, but 
we are able to provide a more incisive condition for the extent of the pairing 
instability in terms of the behavior of the characteristic determinant associ- 
ated with the integral equation (20) (or alternatively, with the characteristic 
determinant of the integral equation (13) for the shape function x{^) [80)117]). 
Concerning the dependence of k^ on effective mass, it is seen in fig. 4 that a 
smaller m* leads to a lower value of the upper critical Fermi momentum (cf. 
further discussion below). 

Numerical results for the OPEC interaction that may be compared with those 
of Takatsuka and Tamagaki and Amundsen and 0stgaard (and with general 
findings of other authors) are displayed in fig. 1 and fig. 5. Fig. 1 refers to 
the uncoupled ^P2 case with m* — 1 and shows results for the kp dependence 
of values of the energy-gap parameter corresponding, respectively, to the 
angle-average approximation (uppermost curve), the pure Af = solution 
(middle curve), and the \M\ = 2 solution (lowest curve). The latter two curves 
are in reasonable agreement with results reported in refs. [57,24,40], having 
maxima near 2 MeV in the vicinity of kp — 2.2 fm~^. It is worth mentioning 
again that our calculations are in accord with the energy ordering of the 
M — and \M\ — 2 solutions found numerically in the early work [24,57,40]; 
moreover, they confirm the generic results that have been derived in sect. 5 
using the separation method and the small parameter expansion. At kp = 
2.2 fm~^ the ratio predicted by eq. (49) is reproduced with better than 1% 
accuracy, and the energy order A^^ < A^^ < A^^ implied by eqs. (49) and (50) 
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is preserved. Considering the near-degeneracies in the condensation energies 
for the different solutions, documented originally by Takatsuka and Tamagaki 
[24], these favorable comparisons warrant some confidence in the application 
of the angle-average approximation to the coupled-channel case. 

Introducing the coupling to the channel, we have solved the triplet pairing 
problem for the OPEG and Argonne potentials in the angle-average approx- 
imation, using the procedure outlined in sect. 5. Fig. 5 shows the calculated 
momentum dependence of the gap components A^^(/c) = Ai(A;) (solid curve) 
and Al^{k) = ^^{k) (dashed curve) for the case of the OPEG potential 
and m* = 1, at a Fermi momentum kp = 2.2 fm~^ near the value yielding 
the maximum strength of triplet pairing. The locations of the maxima, min- 
ima, and nodes of these functions are in close agreement with results given in 
refs. [30,57], and the magnitudes are in general agreement. (It is to be noted 
that in his coupled-channel calculation, Takatsuka included only the M — 
magnetic substate.) 

Figs. 6, 7, and 8 depict the principal results of our coupled-channel •^P2-"^F2 
gap calculations for the more realistic potential models, Vis and v^. Fig. 6 
illustrates the momentum dependence of the L — 1 and L — 3 components 
A|^(/c) for the Argonne vis interaction (abbreviated as Ai(/c) and A2(/c)). As 
expected, the shapes of these gap functions are rather similar to those found 
for the OPEG interaction at the same density (cf. fig. 5). However, the overall 
scale of the pairing effect in the OPEG case is almost 50% larger than for the 
Vis model. 

In fig. 7, we view the dependence of the gap component values A2(A;f) and 
A3(/cf) on the Fermi momentum A-/,', as calculated in angle-average approx- 
imation for the fi8 interaction. Contrasting the results for Ai^kp) with and 
without coupling to the F wave, it is seen that the tensor force acts to sub- 
stantially increase the pairing effect, as has been demonstrated in greater or 
lesser measure by earlier investigators [30,57,40,69]. 

Fig. 8 compares the kp dependence of the ^P2-'^F2 energy-gap parameter Ap 
for the three potentials investigated, the angle-average approximation being 
invoked in all cases. In all three cases, the gap measure A^ peaks in the range 
between kp = 2 and 2.5 fm^^. With the free single-particle spectrum, i.e., 
m* = 1, the maximum value reached by this quantity for the Vu interaction 
is twice that attained for the more refined fig model. As is clear from its 
definition, Ap contains direct contributions both from P-wave and F-wave 
pairing. Inspecting fig. 2, we observe that the greater pairing effect for the 
Vi4 potential relative to Vig can be traced to the greater attraction present 
for Vi4 in both P- and F-wave J = 2 channels over the relevant range of 
the momentum variable k. (The coupling matrix elements Vis{kk) for the 
two potentials are essentially equivalent over the latter range.) Fig. 8 also 
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demonstrates the greater pairing effect of tlic OPEG potential relative to t'lg. 
Such discrepancies between the predictions obtained for different interaction 
models are not unexpected: due to the log singularity intrinsic to the BCS 
pairing phenomenon, the overall strength Ap of the effect in the ^P2-"^F2 
channel is extremely sensitive to the behavior of the input pairing matrix 
elements at momenta in the vicinity of the Fermi surface. In turn, since the 
pairing matrix elements are constructed directly from the given free-space 
NN potential model and since kp ~ (2 — 3) fm~^, one must expect a sensitive 
dependence on properties of the bare interaction that are not determined, or 
only poorly determined, by fits of the NN scattering data up to 350 MeV. 
This point has been developed thoroughly and convincingly in the recent work 
of Baldo et al. [79]. 

In table 1, we report some explicit numerical values of A^^ obtained in our 
calculations based on the two "reahstic" potential models, Argonne vu and 
Vis- A comparison of our numerical results for these potentials with the corre- 
sponding results of Baldo et al. [79] is in order, restricting attention perforce to 
the case of free single-particle energies s{k). First, we point out that the pair- 
ing matrix elements entering our formulas are larger than those appearing in 
ref. [79] by a numerical factor 7r/2, due to different choices of normalization for 
the spherical Bessel functions. Upon accounting for this (immaterial) factor, 
the diagonal matrix elements plotted for the Vis interaction in fig. 3 of ref. [79] 
are apparently in excellent agreement with their counterparts in fig. 2. Second, 
it must be recognized that the quantity D'^{k) defined in eq. (6) of ref. [79] is 
actually twice the function D^{k) appearing in our treatment (see our eqs. (6) 
and (34)). Thus the gap measure Ap = [D^ikp]]^^"^ we have calculated is 
smaller than the Ap quantity of Baldo et al. by a factor 1/V2. The defini- 
tion we have adopted is motivated by the desideratum of consistency with the 
earlier literature [57,40]. In any case, the two calculations should agree once 
this scale factor is applied, since both calculations employ the angle-average 
approximation. 

The results are generally in good agreement; however, there are some modest 
but significant numerical discrepancies. For the ^14 interaction, the current 
results for A^^ versus kp are somewhat higher than those of ref. [79] below 
the peak, and somewhat lower above. The position of the peak, occurring 
near kp = 2.3 fm~^ in the calculation of Baldo et al., is shifted to a slightly 
smaller kp value, near 2.2 fm~^. The value of A^ at the peak is higher in our 
calculation by roughly 5%. In the case of the Vis interaction, the numerical 
agreement is excellent below the peak, but beyond the maximum the current 
results lie below those of ref. [79] by an amount that can reach about 0.1 MeV 
(in the A^;' measure used by Baldo et al.). 

The origin of these relatively minor numerical disparities is not clear. We 
would, however, like to reiterate the concern expressed in refs. [69,79] about 
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the sensitivity of the numerical integration over k in the triplet gap equations 
to the choice of mesh points near the Fermi momentum. The difficulty arises 
because the functions k'^A'l^{k)/E{k) and k'^Al'^{k)/E{k) are very strongly 
peaked at kp- To handle this problem within the present treatment, we have in- 
tegrated out the troublesome log-like part by hand and performed a numerical 
integration only for the remainder, which is much smoother and well-behaved. 
For the latter integration, we employ 288 points, dividing the k domain into 
three intervals: [0, /ci?'], [kF,2kF], and [kpjkjaax], with /cmax = 50 fm~^. Each 
interval contains 96 mesh points distributed in Gaussian fashion - very dense 
around the ends of the interval and sparse in the middle. 

In fig. 8 and table 1, we have included Vis results not only for the baseline 
choice m* = 1, but also for a value m* = 0.78 of the effective- mass parameter 
conventionally thought [30,57] to provide a sensible representation of the in- 
medium single-particle energy. This brings us to the last point that we would 
like to address in the current treatment of triplet pairing in neutron matter. 
Relative to the bare-mass case, the propagator renormalization implied by 
m* = 0.78 yields a substantially smaller gap as well as a lower value for the 
upper critical density (see fig. 4). These effects can be readily understood in 
terms of a reduction of the density of states A'"(0) near the Fermi surface, which 
is related to the effective mass by the familiar expression A^(0) = /cfM^/tt^^^. 
In the singlet-S-wave problem, the approximate formula (67) derived in the 
appendix shows clearly the inverse exponential dependence of the leading com- 
ponent of the gap on the density of states in the case of a negative pairing 
interaction on the Fermi surface, as does the BCS weak-coupling formula in a 
more limited context. The same physical reasoning and similar mathematical 
connections apply in the triplet case. 

It is a common belief that at the higher densities where triplet pairing be- 
comes important, the interactions with particles in the background neutron 
sea produce a considerably greater reduction from the bare mass than is the 
case at inner-crust densities where singlet pairing dominates [57,26]. This con- 
clusion rests on the results of the Brueckner-Hartree-Fock and correlated-basis 
calculations. However, in these calculations, the contribution of pion-exchange 
effects to the renormalization of the nucleon effective mass is ignored. The role 
of the pion degrees of freedom grows rapidly with increasing density, eventually 
giving rise to a rearrangement of the ground state (associated, in particular, 
with pion condensation [118,119]). The latest variational prediction [120] of 
the density Pc at which this phase transition sets in is based on a nuclear 
Hamiltonian containing the Argonne vis two-nucleon interaction (with boost 
corrections to account for the leading relativistic effects) and phenomenologi- 
cal three-nucleon interactions. The result pc ~ 0.2 fm^^ actually lies below the 
density at which triplet pairing is expected to reach its maximum strength. 
Straightforward considerations [121,122] demonstrate that in the vicinity of 
this transition point the contribution of the pion-exchange diagram to the ef- 



29 



fective mass M* diverges and m* goes to infinity. Therefore it is not unfikely 
that the operative value of the effective mass in the domain of triplet pair- 
ing considerably exceeds the bare mass. From this standpoint, conventional 
treatments that assume an m* value somewhat below unity may lose their 
relevance. 



7 Prospectus 

Exploiting the separation approach and the existence of a small parameter, we 
have been able to clarify important general properties of triplet pairing in neu- 
ron matter and to establish a framework and explicit procedures for reliable 
numerical solution of the relevant BCS gap equations. This effort complements 
a recent exhaustive analysis of the full set of generic solutions of the pure (un- 
coupled) ■^P2 problem, also performed with the aid of the separation strategy 
[75]. In view of the corresponding advances that have been made in the ^Sq 
pairing problem [67,80,78], it can be reasonably claimed that the mathemati- 
cal and analytical aspects of BCS pairing of identical fermions in ^Sq and '^P2 
states are now well understood, both within the relevant nucleonic contexts 
and in a more universal, context-independent sense. The same cannot be said 
of the specific physical and mechanistic aspects of triplet or singlet pairing 
of strongly-interacting fermions, which are associated with the realistic in- 
medium particle-particle interaction and renormalized single-particle energies 
that should provide the raw material for the gap equations. A complete and re- 
alistic treatment of pairing in a given strongly-coupled Fermi system demands 
ab initio calculation of these essential inputs. 

The success of the method of correlated-basis (or CBF) theory across a broad 
range of many-body systems and phenomena [124,55] suggests the following 
approach to the quantitative physics of pairing in extended nucleonic systems: 

(a) Dressing of the pairing interaction by Jastrow correlations within CBF 
theory [34,35] 

(b) Dressing of the pairing interaction by dynamical collective effects within 
CBF theory [35,41] (including polarization effects arising from exchange 
of density and spin-density fluctuations, etc.) 

(c) Consistent renormalization of single-particle energies by short- and long- 
range correlations within CBF theory (cf. ref. [123]) 

This approach has already been explored in the ^Sq neutron pairing problem 

[41,55], although the assumed Jastrow correlations have not been optimized 
and only a second-order CBF perturbation treatment is available for step 
(b). Application of this scheme to ^P2-'^F2 pairing in neutron-star matter is a 
challenging but potentially rewarding prospect. 
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The inclusion of other physical effects will also be important in achieving a 
realistic description of pairing in the neutron-star medium. Modifications due 
to three-nucleon interactions should be considered, the calculations should be 
performed in /3-stable matter (cf. ref. [69]), and relativistic effects should be 
examined (see, for example, ref. [70,77]). As pointed out in sect. 6, the impact 
on pairing of mesonic degrees of freedom, especially through an enhancement 
of the nucleon effective mass in the vicinity of pion condensation, deserves 
close investigation. 

In conclusion, it is important to emphasize that although the specific calcu- 
lations performed in the present work were based on pairing matrix elements 
computed directly from bare, in-vacuum NN interactions and on quadratic 
single-particle energies, many of our findings are unaffected by these simplifica- 
tions. Indeed, the analyses that have been performed (leading to the universal 
properties discussed here and in ref. [75]), as well as the calculational proce- 
dures that have been developed, will retain their applicability and essential 
validity when realistic in-medium inputs are employed in the gap equations. 

Within this in mind, further studies in the same spirit - aimed at develop- 
ing accurate computational procedures and uncovering universal properties 
- are of considerable interest. As a continuation of the work of sect. 6, the 
separation-transformed BCS gap equations for pure neutron matter may be 
solved in the ^P2-^F2 coupled channel without resorting to the angle-average 
approximation, but still within the setting of the rapidly convergent small- 
parameter expansion, and with bare pairing matrix elements and unrcnor- 
malized single-particle energies. Another project that awaits conclusion is a 
perturbative treatment of the lifting of the universalities and degeneracies 
that were established for the uncoupled pairing by the analysis carried 
out in ref. [75]. Also, the separation analysis performed in that paper may 
be extended to explore the superfluid phase diagram of liquid "'He - although 
this task is considerably more demanding than in the case of ''P2 pairing in 
neutron matter. Finally, application of the separation approach to studies of 
^Si-'^Di pairing in symmetrical nuclear matter and in the expanding phase of 
heavy-ion collisions may produce new insights into these intriguing problems. 
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9 Appendix 



In this appendix we review the apphcation of the separation method to ^Sq 
pairing [67] and display the first two approximants to the gap amphtude 
Ap = A{kF) within a small-parameter expansion of the separated prob- 
lem [80] . Having set the pattern, we then sketch the extension of the small- 
parameter expansion to pairing. 

To provide a baseline for the separation method, it is useful to revisit the 
problem of a separable pairing matrix V{k, k') = VQ(f){k)(f){k') with (plkp) = 1 
and vq = V'^^^{kF, kp) 7^ 0. The standard gap equation in the ^Sq channel, 



with dro = k dk/27i , becomes 



2^e{k') + A^ik') 



(One may note that the differential volume element dro entering eqs. (55)-(56) 
and other formulas below differs from the volume element dr = {2/7r)k'^dk 
occurring in integrals elsewhere in this paper. The difference arises from a 
different choice of normalization in forming the pairing matrix elements. In this 
appendix, the S-wave pairing matrix elements are defined as in refs. [67,80], 
whereas the convention employed by Takatsuka and Tamagaki [57] for general 
pairing matrix elements is adopted in the rest of the paper.) 

Since the integral on the r.h.s. of eq. (56) is just a number, we see that (/){k) de- 
termines the shape of the gap function, while the integral determines its scale. 
Introducing a dimensionless shape function xik) such that A{k) = Apxik), 
the gap equation for this problem is equivalent to 



X{k) = m, (57) 

l = -vof I '^'^^^ dro . (58) 

J 2^e{k) + Al4>\k) 

Thus, one has only to solve a nonlinear equation for the gap amplitude rather 
a nonlinear singular integral equation for the gap function. We note that a 
solution of this equation exists only when v^? is negative, in which case the 
monotonicity of the r.h.s. implies that there exists one and only one solution. 
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Now we turn to the general situation where V{k, k') need not be separable. 
To gain some benefits of the separable case, let us split the pairing interaction 
into a separable part and a remainder that vanishes when either momentum 
variable is on the Fermi surface: 

V{k, k') = vo(t){k)(t){k') + W{k, k') , (59) 



where W{kF,k) = W{k,kF) = and, as before, (f){k) = V{k. kF)/vF- Again 
writing A(A:) = Apxi^), the gap equation (55) is readily seen to be equivalent 
to the integral equation for the shape function x(/c) 

, f W{k,k')x{k') ^ 



together with the 'algebraic' equation 

1 + Vq — , dro = (61) 

J 2Jeik) + Alx'ik) 



for the gap amplitude Ap (assumed nonzero). It is important to note that since 
W{k, kp) is zero by construction, the integral equation (60) has a nonsingular 
kernel, the log-singularity of the original BCS formulation having been isolated 
in the amplitude equation (61). An iterative solution of this set of equations 
converges very rapidly. Indeed, for neutron matter at densities corresponding 
to the inner-crust region of a neutron star, a good approximation to the gap is 
already obtained by replacing Apxik) in the denominator of (60) by a small 
constant scale factor 5 (e.g. O.OIsf) solving (60) for x{k) by matrix inversion, 
and substituting the result into (61) to obtain the scale factor Ap hy (say) 
Newton's method. 

The existence of a small parameter dp — Ap/sp for the problem plays a key 
role in the success of the separation method based on eqs. (60)-(61). Consider 
the equation (60) for the shape function x{k)- Since the Taylor expansion of 
the remainder function W{k^ k') around k' = kp starts with the first-order 
term dW{k, k')/dk'\k'=kp{k' — kp), the integral over k' in the symmetric vicin- 
ity [kp — 5k, kp + 5k] of this point vanishes for sufficiently smooth functions 
W{k,k') and x{k), 



d'^W{k,k') 



dk' 



dW{k,k') 
5k ^ ^ ' ' 



dk' 



and '^^^5k^l. (62) 
ak 



Also, if beyond this interval we have Apx{k^) -C i{k), or in other words if 5k 3> 
dpkp, it is a safe approximation to replace the A^? term in the denominator 



33 



of (60) by a rather arbitrary scale factor 5. We thereby obtain a practically 
equivalent linear integral equation 

,,,, , f WaJ:')x(k') ^ 



which is always soluble for well-behaved interactions. For most purposes, it is 
permissible simply to set 5 equal to zero. The nature of this "5 approximation," 
and more generally of expansions in the small parameter dp, is developed more 
fully below. 

For realistic nuclear potentials, the accuracy of this approximation in the 
singlet-S problem is in fact much higher than indicated by the simple error 
estimate rfp Inl/dp that accompanies the above preceding. The reason for this 
heightened accuracy is that the main contribution to the integral in (60) comes 
from a region well separated from the Fermi surface. The ^Sq pairing matrix 
elements of realistic nuclear potentials show maxima at ~ 3 fm~^. Hence 
the contribution from the Fermi surface region, already small because of the 
structure of W{k,k'), becomes quite negligible in comparison. 

Based on Taylor-series expansion of the energy denominator of eq. (61), the 
gap amplitude may be systematically expressed as a combination of series 
of terms in {d'j^)"'\nl/dp and (rf|-)" with integral n. Introducing the dimen- 
sionless variable e{k/kpY — 1, eq. (61) is first recast as 



^ . ^Je + 10(e)x(e)de ^ T Ve + 10(e)x(e) ^ ^ ^^^^ 



The integration over e has been divided two intervals, namely (i) the vicinity 
of Fermi surface e G [—1, 1], where terms of both types indicated above will 
appear, and (ii) the tail e e [1, oo], where a simple series in dp is generated. 



J ^/e+l(j){e)x{e)de ^ 7 Ve + 10(e)x(e) 4^2 ? + 10(e)x^(e) ^^ 
/ 2^e2 + 4x'(e) "/ 2|e| ' 2^^ J 2\e\^ 

Extracting the logarithmic contribution from the second term of equation (64), 
we have 



2 , 1 



+ ln^ + ln^ ( 1 + a/1 + 4 



N(0)vf dp 2 

^ij yFTTo(r)A(r) - 1 ^ 1 ?^^^ xATTo(r)x(6) 
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(66) 



The leading approximation to 6?^, 

1 



d% = 2 exp 



,(67) 



is obtained by setting dp = everywhere except within the first logarithmic 
term of eq. (66). The first term inside the square brackets is the usual BCS 
contribution and is clearly dominant if the interaction is weak and localized 
in momentum space around the Fermi surface. If - as is the case for nuclear 
forces - the localization requirement is not fulfilled, the contribution from the 
second term of eq. (67) must also be taken into account. One can then derive 
meaningful results even for potentials with positive vq. The approximant (67) 
is offered as a useful generalization of the BCS weak-coupling formula to a 
broad class of pairing interactions. 

The first correction and higher corrections to the formula (67) may be de- 
veloped as follows. Near the Fermi surface, the denominators appearing in 
eq. (66) are simplified by noting that 



+ 4x'(e) = + 4 + 4 x'(e) - 1 



(68) 



can be expanded in powers of d% [x^(e) — 1] / {e^+d^,) . This is a small quantity: 
for e > dp it is of order dj^/e^, while for e < we can make the expansion 



(69) 



where the subscript stands for e = 0, and use 

2edF 



e2 + d2 - 



< 1 



(70) 



(The prime indicates the derivative with respect to e, while the subscript 
stands for e = 0.) Accordingly, we invoke the expansion 



dl 



1 



ix' - if 

(e^ + dlf 



.(71) 



At this point, any smooth functions in the numerators of the integral terms 
in eq. (66) may be expanded in Taylor series. Since (i) the integrations in the 
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fourth and sixth addends on the l.h.s. of (66) are performed on a symmetric 
interval and (ii) the denominators contain only even powers of e, we need only 
consider terms with even powers of e in the corresponding numerators. For 
example, 



I 



de 



+ dl^e' + dlx\e 



(72) 



\dl 



de- 



(e2 + dp2 



2 ^ 



1 \^ ^ u,p 
I \2 



ix'oT + Xo 



-2 + 2 In 



de 



df 



[t' + dl)-- 
+ ?44(x'o 



/ \2 



2 1 
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This step reveals a slight complication in the analysis. The integral of the first 
term of the Taylor expansion of — 1 vanishes, and consequently the integral 
of the second term of this Taylor expansion is of the same order as the integral 
of the term of the expansion (71) containing [y^ — 1)^ are of the same order 
in dp- 

In evaluating corrections to the result (67), it is convenient to define 



1 



^2n 



/'(6^+4r^ 



(73) 



and 



A„=(x^TT«- 



(74) 



The first correction to d^p is then determined by applying recursion relations 
involving the (see ref. [117]). Eq. (64) is written as 



d% " 



1 

4' 



^4^) In 



2 

dp ' 



(75) 



where 



B 



oc 

I 



0(6)x'(6)de 



-1 



1 



V^0X' - 1 - ^e'As 



de 
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--A, + -As + -A, + 1, 



(76) 



and the improved estimate is 



1 + 7(4) [B-Asln 



■F 



) 



(77) 



follows. 



A calculation for the Argonne V14 interaction illustrates the utility and accu- 
racy of the leading approximation (67) to the energy gap A^^ = SpdF within 
the small-parameter expansion. In fig. (9), the kp dependence of as given 
by the zeroth-order approximation in dp is compared with the essentially exact 
result obtained by iteration of eqs. (60) and (61). Free single-particle energies 
(m* = 1) are assumed. The estimate (67) produces excellent results, especially 
at densities beyond the peak of the A^r curve. 

The same strategy of constructing approximants to the gap problem through 
small-parameter expansion is applicable to pairing in higher angular momen- 
tum states and in particular to pairing. By way of illustration, consider 
the Type 1 (M = ±2) solution of the uncoupled-channel problem. The trans- 
formed gap equation (28) associated with this case has a form similar to that 
for ^Sq pairing, with differences arising from the anisotropy of the P-wave 
problem. A weight factor 3(sin^^)/2 = 3(1 — cos^^^)/2 appears in the inte- 
gration over angles in eq. (28), and 3A|,(1 — cos^ 9)/2 appears in place of A|. 
inside the square-root energy denominator. Doing the k integral first, we can 
directly exploit results and relations of sects. 4 and 5 (specifically, the first of 
eqs. (25) along with eqs. (26) and (36)). Ignoring corrections proportional to 
d'^p = Ap/ep, the only additional term that needs to be considered, relative to 
the angle-average approximation, is 



For the Type 2 solution (with M — 0), the additional term has the shghtly 
different expression (see the last of eqs. (25)) 



The most salient points are that both these contributions are of order unity 

relative to the dominant logarithm in the gap and both are negative. We may 
thus conclude that the true gap will be smaller than that given by the angle- 
average approximation and that the correction to the gap value calculated in 
this approximation will generally be small. 




(78) 




(79) 
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Table 1 

Energy gap Ap for ^P2-^F2 pairing in neutron matter at several values of the Fermi 
momentum, for the Argonne vi^ and V14 interactions and the OPEG potential. 
Free single-particle energies (m* = 1) are used, except for the last column, where 
m* = 0.78. 



kp 


OPEG 


''11 


''LS 


ri,s(0.78) 


1.0 


0.0009 


0.0020 


0.0020 




1.5 


0.1369 


0.2496 


0.1576 


0.0263 


2.0 


0.5376 


0.8028 


0.4362 


0.0947 


2.2 


0.6167 


0.8519 


0.4235 


0.0875 


2.3 


0.6175 


0.8182 


0.3807 


0.0740 


2.5 


0.5419 


0.6542 


0.2458 


0.0384 


3.0 


0.1538 


0.1330 


0.0059 





0.4 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Fig. 1. Test of the angle-average approximation for ^P2 pairing, without coupling 

to the ^F2 state. The energy gap parameter Ap (square root of the quantity (7), 
evaluated at the Fermi surface) is plotted versus Fermi momentum kp, for the pure 
M = solution (long-dashed curve) and the |M| = 2 solution (short-dashed curve), 
and as calculated in the angle-average approximation (solid curve). The pairing 
matrix elements are constructed from the OPEG potential, and free single-particle 
energies are assumed (m* = 1). 
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Fig. 2. Diagonal pairing matrix elements Vl^i^{k,k) = V\i{kk) (solid curves), 
Viz{kk) (dashed curves), and Vzz{kk) (dot-dashed curves) in the ^P2-^F2 coupled 
channel, as functions of momentum k. Curves for the Argonne vn and uig interac- 
tions are labeled with a square and with a diamond, respectively. 
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Fig. 3. Illustration of gap existence and closure for pairing by the Argonne 
Vis interaction, with coupling to the state turned off. Dependences on Fermi 
momentum kp, of the diagonal matrix elements Vi{^'^{kF, kp) = Vup of the pair- 
ing interaction (solid curve), the dimcnsionless characteristic determinant of the 
linearized version of system (20) (dashed curve), and the gap component value 
^i^ikp) = determined in angle-average approximation (dotted curve with 
diamonds). Free single-particle energies are assumed. 
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Fig. 4. The characteristic determinant for pairing in the ^P2-^F2 coupled channel as 

a function of Fermi momentum kp. Results are shown for the Argonnc vis potential 
and three values of the parameter m* = M*/Mn in an effective-mass approximation 
for the normal-state single-particle energies: m* = 1 (short-dashed curve), m* = 0.8 
(long-dashed curve), and m* = 0.6 (solid curve). 
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Fig. 5. Gap components A'l^^ (k) = Ai{k) (solid curve) and A|^(A;) = A3(/c) 
(dashed curve) for ^P2-'^F2 pairing (coupled-channel case), calculated in an- 
gle-average approximation for the OPEG potential and m* = 1, and plotted versus 
momentum k at kp = 2.2 fm~^. 
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Fig. 6. Gap components Ap^(A;) = Ai{k) (solid curve) and A|^(A;) = A3(/c) 
(dashed curve) for ^P2-'^F2 pairing (coupled-channel case), calculated in an- 
gle-average approximation for the Argonne vis potential and m* = 1, and plotted 
versus momentum k at kp = '2.2 fm~^. 
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Fig. 7. Gap-component values Al^{kF) = Ai(fci^) (solid curve) and 

A'l^{kp) = A-i{kp) (long-dashcd curve) with channel coupling, and Ai(A;f) without 
channel coupling (short-dashed curve), calculated in angle- average approximation 
for the Argonne vis potential and m* = 1, and plotted as a function of Fermi 
momentum kp. 
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Fig. 8. Dependence on Fermi momentum of the overall gap measure Ap (square 
root of the quantity (7), evaluated at the Fermi surface) for coupled-channel ^P2-^F2 
neutron pairing, as determined in the angle-average approximation. Dotted curve: 
for Argonne ?;i4 interaction and effective- mass parameter m* = 1. Short-dashed 
curve: for OPEG potential and m* = 1. Solid curve: for Argonne vis potential and 
m* = 1. Long-dashed curve: for Argonne vig interaction and m* = 0.78. 
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Fig. 9. Superfluid energy-gap amplitude Aj? in the So channel, plotted as a function 
of Fermi momentum kp for the Argonne vis potential and m* = 1, as calculated 
"exactly" by iteration (solid curve) and to leading approximation in the smallness 
parameter dp (dashed curve). 
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